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ABSTRACT 

We study the vertical gradient in azimuthal velocity of spiral galaxy NGC 4244 in 
a thin disk model. With surface density accounting for the rotation curve, we model 
the gradient properties in the approximation of quasi-circular orbits and find the 
predictions to be consistent with the gradient properties inferred from measurements. 
This consistency may suggest that the mass distribution in this galaxy is flattened. 
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1 INTRODUCTION 


Measurements of galactic kinematics carries the information 
of how gravitating mass is distributed in galactic interiors. A 
basic quantity encoding this information is a rotation curve. 
But for a spiral galaxy appropriately aligned with respect 
to the observer, it is also possible to ascertain and study the 
vertical structure of the rotation above the mid-plane which 
provides a further piece of information indispensable for a 
more precise determination of the galactic mass distribution. 

Measurements enabling a more detailed detection of 
the velocity field structure have be en performed, among 
others, for spiral galaxies: NGC 891 (lOoste rloo et al.ll2007[ 


i pirai galaxies: INUC biJl (lUoste rloo et m.l l2UU<1 : 
I2OO6II NG C 4559 tearbieri et al.ll2005ll. N GC 


iHeald et al. _ _ _ _ _ _ .. _ _ _ 

4302 jHeald et al.l^2007^. NGC 577 5 (|Heald et al.ll200fill . for 


NGC 2 403 llFr aternali et al.l [2002fl and for the Milky Way 
galaxy llLevine et al.l 2008h . In particular, it has been estab¬ 
lished that the azimuthal component of the rotation gradu¬ 
ally diminishes in almost a linear fashion with the altitude 
above the mid-plane. 

Recently, we tested the hypothesis that spiral galaxies, 
in contrast to what seems to be true for most galaxies, do not 
necessarily have to be dominated by massive non-baryonic 
dark matter halos - at least it cannot be excluded that some 
of these galaxies might be flattened disk-like objects. But, in 
the absence of dominating spherical components, the motion 
of matter in the galactic mid-plane vicinity, to within the ap¬ 
proximation of axial symmetry, should be satisfactorily well 
described as if governed entirely by the gravitational field of 
a thin (or a finite-width) axi-symmetric disk with surface (or 
column) density accounting for the observed fragment of the 
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rotation curve. From a good mass model it is additionally 
required that its predictions should be consistent with other 
observations, in particular, the model should predict correct 
behavior in the structure of the galactic rotation held. In 
this context we may use the vertical gradient observable as 
a test of the disk model of spiral galaxies. 

We dealt with the v ertical gradient in d isk model in 
lljalocha et al.|[2010lj and lljalocha et al.|[201lli . To describe 
the gradient in this model we additionally assumed that the 
projections of orbits onto the mid-plane were approximately 
circular, which is the subject of the so called quasi-circular 
orbits approximation assumed also in the present paper. In 
lljalocha et al.ll2010lj we also verihed our approach by per¬ 
forming a numerical simulation of motion of test bodies on 
various orbits in the gravitational potential of a thin disk and 
we compared the analytical estimates of our approach with 
the simulated fall-off in azimuthal velocity near the mid¬ 
plane. The results of the simulation and of our approach were 
consistent with each other. In the thin disk model framework 
we were able to reproduce the observed values of the gradi¬ 
ent and its properties, such as a linear decrease of azimuthal 
velocity with the altitude above the mid-plane. We consider 
this as a success of the thin disk model. 

Furthermore, in lljalocha et al.l[^014l l we applied a more 
general, finite-width disk model to study the vertical gradi¬ 
ent in the Milky Way. The width of the disk was constrained 
so as to obtain the best fit to microlensing measurements 
being another source of information about the mass dis¬ 
tribution, quite independent of the rotation curve. On this 
occasion we had the opportunity to test limits on the ap¬ 
plicability of the thin disk to modeling the vertical gradient 
of rotation. The conclusion was intuitively clear: when the 
gradient measurements were performed not too close to the 
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mid-plane, outside the main concentration of mass of the 
finite-width disk (for heights at least 0.5 kpc in the case of 
the Milky Way galaxy), then the approximation of the in¬ 
finitesimally thin disk was sufficient and there was no need 
for using the computationally more demanding finite-width 
disk. 

The two observations above are the basis for our use of 
the thin disk modeling of the vertical gradient in the present 
paper. We expect that the approach could be applicable also 
to other flattened galaxies. 

NGC 4244 is another galaxy for which the verti- 
cal structure in the rot ational velocity can be studied. 
IZschaechner et ^ ll201lll published a velocity field above 
the mid-plane as well as the resulting rotation curve in both 
galactic halves. The reported results for the detected gradi¬ 
ent in their numerical model of NGC 4244 were comparable 
in both galactic halves: —9^2 approaching half 

and —9'^2 receding half, with the gradient de¬ 
creasing in magnitude to —5^2 ~^-2 s.*kpc ^ 

radius of 10 kpc in these halves, respectively. The authors 
pointed out that the lag in Hi was present in absence of an 
Hi halo. This is an indication that the mass distribution in 
NGC 4244 may be flattened, and is what motivates our ap¬ 
plying of the thin disk to model this particular galaxy. We 
want to check whether our simple model applied to NGC 
4244 would account for the reported gradient values and 
their variability with the radial distance. 


2 MASS DISTRIBUTION IN THIN DISK 

In disk model we are looking for a surface density accounting 
for the entire observed rotation curve under the assumption 
of axial and mid-plane-reflection symmetries. 

For obvious reasons, rotation curves do not extend far 
enough and one has to extrapolate the rotation data. But 
then, the model results unavoidably depend on the way one 
chooses to extrapolate beyond the last measured point. This 
feature is not a failure of a mass model as such but is inherent 
in the gravitation of flattened mass distributions. In general, 
for a flattened mass distribution (if not representable as a 
union of contiguous similar homeoids), the rotation value at 
a given radius is depended also on the masses distributed 
exterior to that radius, likewise, column mass density at a 
given point is a point-dependent integral functional of the 
entire rotation curve. Therefore, the thin disk, as a model of 
a flattened mass distribution, must be used with due care. 
A more detailed discu ssion of these issues were described 
by iBratek et al.l (|2008l 'l. Sometimes, it is possible to find a 
global solution by iterations when additional measurements 
complementary to the rotation d ata are known beyon d the 
last measurement of the rotation (|,Talocha et ahll^OOSh . 

2.1 Rotation curve of NGC 4244 

As for NGC 4244, we take the rotation curve as the start¬ 
ing point. We make use of t h e rota tion measurements pub¬ 
lished by IZschaechner et ^ (l201ll ') . They were carried out 
separately for both halves of the galaxy, but we have super¬ 
imposed and averaged out these data because the modeling 
should be done consistently with the axial symmetry the 


disk model approximation postulates. This approximation 
should suffice to model the vertical gradient to the approx¬ 
imation of axial symmetry, because the gradient values in¬ 
ferred from measurements, as it was mentioned earlier, are 
comparable in both halves. 

The rotation curve which we adopt in this work is shown 
in figure Fig[T] It is obtained by interpolating the measured 
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Figure 1. Circular velocity of NGC 4244 in the mid-plane as 
a function of the galactocentric distance (obtained by assuming 
a distance to NGC 4244 of 4.4 Mpc). The collection of points 
show the rotation values in the approaching [black circles] and 
receding [open circles] halves read from curves published in 
llZschaechner et al ][2oTi . The [solid line] is our adopted rotation 
curve obtained by means of interpolating between the rotation 
values in both halves and extended beyond ^ 17 kpc by adjoining 
a Keplerian tail. 

rotation values and, additionally, by extending them beyond 
« 17 kpc by a Keplerian tail. The purpose of this artificial 
extension is only to facilitate integration to infinity in find¬ 
ing the corresponding surface density, and it will not signif¬ 
icantly affect the surface density at lower radii in the region 
of interest for the gradient estimation. This extension is not 
entirely arbitrary, however. A decrease in a Keplerian fash¬ 
ion in the rotation, in a neighborhood of the last measured 
point, was re ported in another work investigating NGG 4244 
ll011inalll996lf . This observation may be used to support the 
form of the artificial extension we adopted. 


2.2 Surface density in NGC 4244 

In modeling the mass distribution and the vertical gradient, 
we limit ourselves to using the thin disk model, instead of 
applying a fin i te-wid th disk, because, as it was illustrated by 
Ijalocha et al.l (l2014l f . when the gradient is calculated above 
the mid-plane, outside the main concentration of masses 
characterized by some vertical width-scale, then the vertical 
structure becomes insignificant for the gradient determina¬ 
tion in this region and one can use a simpler model. 

We determine a surface density in thin disk model based 
on a given circular velocity, by using a formula turning 
the velocity squared to the surface density as derived by 
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iBratek et alJ (l2014l l: 

( 1 ) 

P 

This integral should be evaluated in the principal value sense 
about the singularity a; = 1 of the expression in the square 
bracketsQ 

As we have already pointed this out, the surface mass 
density is uncertain close to the radial position of the last 
point on the rotation curve and beyond it, owing to to the 
nature of a flattened mass distribution as such and lack 
of measurement data. In some cases we could overcome 
this indeterminacy by additionally using a surface density 
of gas distributed beyond the last point and find a global 
rot ation by itera t ions in a way similar to that proposed 
by Ijalocha et al.l ll2008h . But there is no such complemen¬ 
tary measurements available for NGC 4244, and we have to 
choose a different approach. 

Since we are interested in the gradient modeling in the 
galactic interior r < lOkpc separated well enough from the 
last measured point at « 17kpc, we chose to extrapolate 
the rotation by smoothly joining to it a Keplerian tail, so 
as to extend the natural slope of the curve to larger radii, 
as shown in Fig[T] The internal part of the corresponding 
surface density is not significantly dependent on this exten¬ 
sion to within reasonable distortions of that extension. To 
support this statem ent, we recall a calculation presented in 
dBratek et al.ll20(^ . which implies that in thin disk model 
the relative contribution to the local surface density from an 
almost Keplerian tail can be neglected at radii lower than 
roughly 2/3 of the radial extent of the measured part of the 
rotation curve, provided the mass function of the flattened 
mass component can be assumed to be almost saturated 
at the last measurement point. Then, in practice, it would 
make no difference if we set the rotation to be zero beyond 
the last measured point or not. 

The surface density corresponding to the adopted rota¬ 
tion curve shown in Fig[T]is presented in Fig[3 It should be 
stressed that in disk model, unlike in spherical symmetry, 
the Keplerian tail in a region is not equivalent to lack of 
matter in that region, which is readily seen in Fig[2] where 
surface density is indeed non-vanishing in the region beyond 
« 17kpc, where the rotation curve was set to be Keplerian. 
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3 VERTICAL GRADIENT IN AZIMUTHAL 
VELOCITY 

In th e quasi-circular orbits approximation djalocha et al.l 
l201Ct) . the azimuthal velocity at an altitude 2 ; above the 


^ We use complete e l liptic integrals K and E in forms as defined 
bv iGradshtein et al.l ll2007h : 


7r/2 

K{k) = J 
0 


\/l — sin^ (j) 


77/2 

E{k) = J 
0 


Note, that they differ in the a dopted convention for the arg ument 
notation from those given bv lAbramowitz Sz StegunI lIlQTj') . 
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Figure 2. Surface mass density accounting in thin disk model 
for the adopted rotation curve of galaxy NGC 4244 shown in 
Fig^ Unlike in spherical symmetry, surface density is nonzero in 
a region r > 17 kpc with Keplerian rotation curve. 


mid-plane is found from 


„5(r,*) = /-/Eskiai 

iyir + xf + 



2 2,2 
X -r +z 

(r - xf + 



( 2 ) 


4rx 


(r -I- x)^ + ’ 


where the surface density expressed by Eq[T]have been sub¬ 
stituted for cr(x). This relation is used to describe the behav¬ 
ior of the azimuthal velocity component in a region above 
the mid-plane and exterior to the main concentration of 
masses of the flattened mass distribution, where the ver¬ 
tical structure of this distribution can be neglected, but still 
close enough to the mid-plane, so that the projection of or¬ 
bits onto the mid-plane could be regarded as a perturbation 
of the concentric circular orbits based on which Eq[T] turn¬ 
ing the circular velocity to the surface density (and its in- 
verse) are derived in thin disk model (though, a reasoning in 
jjalocha et al.ll2010ll admits a more general class of orbits). 


We will estimate the vertical gradient in azimuthal ve¬ 
locity with the help of two methods. The first method gives a 
mean magnitude of the gradient by means of a linear approx¬ 
imation of the azimuthal velocity. The mean value, denoted 
by 7 , is determined in a given rectangular region (r, z) by 
finding a best fitting profile of the form 


v^{r, z) = v^{r,0) + X \z\ + Sz, (3) 


approximating the values of azimuthal velocity stipulated at 
several points arrange d in an array i n that region. The linear 
formula was used by iLevine et al.l (l2008l i to estimate the 
vertical fall-off rate in the rotation speed of the Milky-Way 
galaxy. Here, the use of the linear expansion is justified by 
the linearity of the decrease in azimuthal velocity predicted 
by the thin disk model for NGC 4244, as illustrated in FigO 
Owing to this behavior (and the reflection symmetry with 
respect to the mid-plane z = 0, which implies that the rolling 
parameter S must be set equal to zero in disk model), it 
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Figure 3. Behavior of azimuthal velocity off the mid-plane in 
galaxy NGC 4244, predicted in quasi-circular orbits approxima¬ 
tion in thin disk model (that is, calculated with the help of inte¬ 
gral Eq[2]l. The thin lines are curves of constant radial variable: 
r = 8 , 9,10,11,12kpc, counting from the bottom line to the top 
line. Each point represents the mean azimuthal velocity compo¬ 
nent, with the respective bar representing the standard deviation. 
The velocity decreases almost linearly with the altitude \z\ above 
the mid-plane. 



5r 


6z 

S-V4, r km 1 
dz ^ s- kpc J 

re(l, 16) 

0.2 

2e(o.5,1) 

0,02 

-8.38 ±0.33 

re (2,10) 

0.2 

2e(0.5,3) 

0.2 

-8.24 ±0.21 

re (10,16) 

0.2 

26(0.5,1) 

0.02 

-5.23 ± 0.05 

rG(10,11) 

0.01 

2 6(0.5,3) 

0.2 

-6.00 ±0.00 


Table 1. Rectangular arrays of rotation data (z, r , Sz and 5r, are 
expressed in kpc) and the corresponding mean vertical gradient 
of rotation (coefficient 7 in EqlSj. 


Fig[4] illustrates the gradient dependence in function of the 



Figure 4. Absolute values of the vertical gradient in azimuthal 
velocity in thin disk model shown in function of the altitude z 
above the mid-plane for various values of the radial variable from 
within the interval r G (2,16) kpc in steps of 2 kpc. These lines 
represent the values obtained using the integral expression Eq[4l 

altitude above the mid-plane for various radii in a region 
from 2 kpc to 16 kpc. The absolute gradient value depends 
weakly on the 2 variable, except close to the galactic center. 
This reflects the linear decrease of azimuthal velocity evident 
in FigO The maximum gradient magnitude exceeds even 
20 but only close to the center and at low heights above 
the mid-plane. In a prevailing part of the galactic region 
the gradient magnitude does not exceed 10 , decreasing 

moderately with the radial distance and is lower than 4 
for the outermost radii (compare, the r — 16 kpc line in 
FigH}. 

FigO shows for several altitudes the variation in the 


is rational to consider only a single number describing the 
vertical velocity fall-off for NGC 4244, namely, the mean 
gradient magnitude 7 . 

Table [T] shows the results of applying the first method 
to the theoretical rotation values predicted by evaluating 
the integral Eq[2] inside several rectangular regions defined 
in that table. To some degree, the gradient results depend 
on the horizontal and vertical extent of a given region where 
the gradient is estimated. To trace this dependence we use 
a second method. 


The second method of determining the vertical gradient 
in azimuthal velocity uses a direct expression for the gradi¬ 
ent, obtained by differentiating the expression for velocity 
Eq[2] with respect to 2 and using some identities between 
elliptic integrals K, E and their derivatives. The result is 


dzV^ (r, 2 ) 


Gz r x<^(x)dx 


(4) 




yJ + {r- x) 


( 2 ^ + {r-x) ) 



Figure 5. Absolute values of the vertical gradient in azimuthal 
velocity shown as functions of the radial variable for various alti¬ 
tudes above the mid-plane. 


gradient in function of the radial distance. In the central 



























5 


part the absolute gradient value is not only high (exceed¬ 
ing 20 for low altitudes) but it also strongly changes 
with the radial variable. But starting from a radius of 4kpc 
the gradient magnitudes are much lower (below 10 
and they exhibit a significantly weaker dependence on the 
radial variable, being more pronounced only for lower alti¬ 
tudes above the mid-plane. 

We can summarize the above results by saying that the 
vertical gradient predicted for NGC 4244 is not too high - 
its mean value in a rectangular region r £ (1,16) kpc and 
\z\ £ (0.5, 1) kpc is —8.38±0.33 , and it would not differ 

much from a value of —8.24 ± 0.21 obtained by shrink- 

ing the radial interval to r £ (2,10) kpc and increasing the 
height above the mid-plane to 3 kpc. We predict also a de¬ 
crease in the gradient magnitude for higher radial distances: 
in a rectangle r £ (2,16) kpc and \z\ £ (0.5, 1) kpc we obtain 
a mean gradient value of —5.23 ± 0.05 ■ These predic¬ 

tions are in agreement to within error limits with t he gradi¬ 
ent properties reported in llZschaechner et ^1201 il l, namely, 
with the gradient value of —ffi® approaching 
half and —in the receding half, and with the gra¬ 
dient magnitude found to be decreasing to 5^2 nnd 

'^-2 near a radius of 10 kpc in the approaching and 

receding halves, respectively. 

In the first method we defined the error of our model 
gradient value as the uncertainty in the parameter 7 inferred 
from a linear regression model Eq[3] fitting to the velocity 
values predicted for quasi-circular orbits. The error does not 
include the error transferred from the uncertainties in the 
rotation curve determination. We can get a good impression 
of the scale of these uncertainties from Fig|^ which illus¬ 
trates the influence of variations in the rotation curve on 
the predicted gradient val ues. But the relative erro rs in the 
gradient determination in llZschaechner et ^l201ll ) are even 
higher, therefore a more precise analysis of the model un¬ 
certainties seems of minor importance for our conclusions 
concerning the gradient properties. 


4 CONCLUDING REMARKS 

There are various models aimed at explaining the presence 
and behavior of the vert i cal gradient of azimuthal ve locity 
(e.g.. lCollins et all (|2002l '): lFraternali fc Binn^ (|2006l 'l'l. but 
our model is characterized by high simplicity and minimum 
of assumptions. We use a disk model with a surface den¬ 
sity accounting for the full rotation curve measured for a 
galaxy, assuming that the motion of matter is entirely due 
to the gravitational potential of the disk. This means that 
most of matter is assumed to form a flattened mass dis¬ 
tribution, without a dominating heavy halo. Despite this 
simplicity our model accounts for the vertical gradient in 
azimuthal velocity, predicting correct gradient magnitudes 
and its approximate independence of the altitude above the 
mid-plane. 

NG C 4244 is the next among t he galaxies examined by 
us so far (j.lalocha et al.l (|201fll . [201ll B for which our approach 
to modeling the gradient turned out to perform satisfacto¬ 
rily well, predicting vertical gradient magnitudes and their 
behavior consistently with expectations based on the mea¬ 
surements. The most controversial in this approach may be 
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Figure 6. This figure shows results analogous to those in Fig[5] 
for various rotation curves. Thick lines correspond to a rotation 
curve obtained by smoothing spline interpolating between the 
approaching- and receding- half rotation curves (it almost over¬ 
laps the adopted rotation curve in Fig[T] which shows that the 
averaging method is not influential to the overall shape of that 
curve). Thin lines correspond to deformed rotation curves ob¬ 
tained by taking linear combinations aVa + (1 — a)vr of the rota¬ 
tion curve in the approaching half (denoted here by Va) and that 
in the receding half (denoted by Vr). 

our not including the other mass components except the 
disk-like only. But obtaining satisfactory results in model¬ 
ing a galaxy in this approach, simply may indicate that the 
galaxy has a dominating flattened mass component, at least 
in the region where the gradient is measured and modeled. 
This we give as the main conclusion of our work. 

An interesting property we want to bring to the at¬ 
tention are high gradient magnitudes in the galactic center 
vicinity of NGG 4244 (mainly close to the mid-plane), ex¬ 
ceeding the mean gradient magnitude determined for this 
galaxy. A similar pheno menon occurs for another galaxy. 
IZschaechner et ahl (gOlJ) mention that the measured gradi¬ 
ent magnitude in the central part of NGG 891 is very high 
(reaching 43 ^:^^^) although it decreases strongly for outer 
radii (to 14 This may suggest that high gradient mag¬ 

nitudes in the central parts of spiral galaxies (much higher 
than the mean magnitude for the whole disk) might be a 
qualitative property of the vertical gradient which the disk 
model accounts for well. 
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